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Abstract. The aim of this work is to enumerate alternating sign matrices (ASM) that are quasi-invariant under a 
quarter-turn. The enumeration formula (conjectured by Duchon) involves, as a product of three terms, the number of 
unrestricted ASM's and the number of half-turn symmetric ASM's. 

Resume. L'objet de ce travail est d'enumerer les matrices a signes alternants (ASM) quasi-invariantes par rotation 
d'un quart-de-tour. La formule d' enumeration, conjecturee par Duchon, fait apparaitre trois facteurs, comprenant le 
nombre d' ASM quelconques et le nombre d'ASM invariantes par demi-tour. 

1 Introduction 

An alternating sign matrix is a square matrix with entries in { — 1, 0, 1} and such that in any row and 
column: the non-zero entries ahernate in sign, and their sum is equal to 1. Their enumeration formula was 
conjectured by Mills, Robbins and Rumsey (||), and proved by Zeilberger (^, and almost simultaneously 
by Kuperberg (^. Kuperberg used a bijection between the ASM's and the states of a statistical square- 
ice model, for which he studied and computed the partition function. He also used this method in to 
obtain many enumeration or equinumeration results for various symmetry classes of ASM's, most of them 
having been conjectured by Robbins (^. Among these results can be found the following remarkable one. 

Theorem 1 (Kuperberg). The number Aq^{AN) of ASM's of size AN invariant under a quarter-turn 
( QTASM's) is related to the number A{N) of( unrestricted) ASM's of size N and to the number Ajj'p (2-/V) 
ofASM's of size 2N invariant under a half-turn (HTASM's) by the formula: 

Aqt(47V) = A^j{2N)A{N)\ (1) 

More recently, Razumov and Stroganov (^ applied Kuperberg's strategy to settle the following result 
relative to QTASM's of odd size, also conjectured by Robbins (0) . 
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Theorem 2 (Razumov, Stroganov). The numbers of QTASM's of odd size are given by the following 
formulas, where Ajj'p(2A^ + 1) is the number ofHTASM's of size 2N + 1: 

Aqj{4N-1) = Ai^j{2N ~ 1)A{N)'^ (2) 
Aqjim+l) = Aiij{2N + 1)A{N)\ (3) 

It is easy to observe (and will be proved in Section^ that the set of QTASM's of size 4N + 2 is empty. 
But, by slightly relaxing the symmetry condition at the center of the matrix, Duchon introduced in (^ 
the notion of ASM's quasi-invariant under a quarter turn (the definition will be given in Section |^ whose 
class is non-empty in size AN + 2. Moreover, he conjectured for these qQTASM's an enumeration formula 
that perfectly completes the three previous enumeration results on QTASM. It is the aim of this paper to 
establish this formula. 

Theorem 3 The number Aqji^N + 2) ofqQTASM of size 4N + 2 is given by: 

Aqj{4:N + 2) = Anj{2N + 1)A{N)A{N + 1). (4) 

This paper is organized as follows: in Section ^, we define qQTASM's; in Section |^, we recall the 
definitions of square ice models, precise the parameters and the partition functions that we shall study, 
and give the formula corresponding to equation (Q) at the level of partition functions; Section ^ is devoted 
to the proofs. 

2 ASM's quasi-invariant under a quarter-turn 

The class of ASM's invariant under a rotation by a quarter-turn (QTASM) is non-empty in size AN — 1, 
AN, and AN + 1. But this is not the case in size AN + 2. 

Lemma 4 There is no QTASM of size AN + 2. 

Proof: Let us suppose that M is a QTASM of even size 2L. Now we use the fact that the size of an ASM 
is given by the sum of its entries, and the symmetry of M to write: 

2L^ ^ Af,,,=4x ^'kj (5) 

l<ij<2L l<ij<L 

which implies that the size of M has to be a multiple of 4. □ 

Duchon introduced in (^ a notion of ASM's quasi-invariant under a quarter-turn, by slightly relaxing 
the symmetry condition at the center of the matrix. The definition is more simple when considering the 
height matrix associated to the ASM, but can also be given directly. 

Definition 5 An ASM M of size AN + 2 is said to be quasi-invariant under a quarter-turn (qQTASM) if its 
entries satisfy the quarter-turn symmetry 

M4]y+2-j+l,4N+2-i+l = Mij (6) 



except for the four central entries {M2n,2N, M2n,2N+i, M2n+i,2N, M2n+i,2N+i) that have to be either 
(0,-1, -l,0)or(l, 0,0,1). 



Quarter-turn Symmetric Alternating Sign Matrices 



51 



We give below two examples of qQTASM's of size 6, with the two possible patterns at the center. 
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In the next section, we associate square ice models to ASM's with various types of symmetry. 

3 Square ice models and partition functions 

3. 1 Notations 

Using Kuperberg's method we introduce square ice models associated to ASM's, HTASM's and (q)QTASM's. 
We recall here the main definitions and refer to (Q) for details and many examples. 

Let a G C be a global parameter. For any complex number x different from zero, we denote x = 1/x, 
and we define: 

(t(x) ~ x — X. (7) 

If G is a tetravalent graph, an ice state of G is an orientation of the edges such that every tetravalent 
vertex has exactly two incoming and two outgoing edges. 

A parameter a; ^ is assigned to any tetravalent vertex of the graph G. Then this vertex gets a weight, 
which depends on its orientations, as shown on Figure ll|. 



a{a'^) o'(a^) (T{ax) a{ax) a{ax) a{ax) 
1-10 



Fig. 1: The 6 possible orientations, their associated weights and the corresponding entries in ASM's 

It is sometimes easier to assign parameters, not to each vertex of the graph, but to the lines that compose 
the graph. In this case, the weight of a vertex is defined as: 



y 



xy 



When this convention is used, a parameter explicitly written at a vertex replaces the quotient of the 
parameters of the lines. 

We will put a dotted line to mean that the parameter of a line is different on the two sides of the dotted 
line. We will also use divalent vertices, and in this case the two edges have to be both in, or both out, and 
the corresponding weight is 1: 



1 1 



52 



Jean-Christophe Aval, Philippe Duchon 



The partition function of a given ice model is then defined as the sum over all its states of the product 
of the weights of the vertices. 

To simplify notations, we will denote by Xjv the vector of variables [xi , . . . ,xn)- We use the notation 
X\x to denote the vector X without the variable x. 

3.2 Partition functions for classes of ASM's 

We give in Figures ^ ||, and ^ the ice models corresponding to the classes of ASM's that we shall study, 
and their partition functions. The bijection between (unrestricted) ASM's and states of the square ice 
model with "domain wall boundary" is now well-known {cf. (Q)), and the bijections for the other symmetry 
classes may be easily checked in the same way. The correspondence between orientations of the ice model 
and entries of ASM's is given in Figure |lj 



XN- 



Z{N;xi, . . . ,xn,xn+i, ■ ■ ■ ,X2n) = ^ 

X2^ 

Xl^ ^ 

XN+1 X2N 

Fig. 2: Partition function for ASM's of size A'^ 

With these notations. Theorem |^ will be a consequence of the following one which addresses the 
concerned partition functions. 

Theorem 6 When a = ujq ^ exp(i7r/3), one has for N > I: 

Zqj{AN; X2N-UX, y) = a(a)-i%x(2^; X2n-ux, y)Z{N; X2N-1, x)Z{N; X2N-i,y) (8) 

and 

Zqj{4N + 2; X2N, X, y) = a{a)-^ Z-^^{2N + 1; X2N. x, y)Z{N; X2n)Z{N + 1; X2N.X, y). (9) 

Equation is new; equation (||) is due to Kuperberg for the case x = y. To see that Theorem 
^implies Theorem |](and Theorem [l|), we just have to observe that when a — uif^ and all the variables 
are set to 1, then the weight at each vertex is cr(a) = cr(a^) thus the partition function reduces (up to 
multiplication by (7(0)""™''°'' °* vorticcs-j number of states. 

4 Proofs 

In this extended abstract, we shall only give the main ideas of the proofs. Most of them are greatly inspired 
from (^. To prove Theorem ^, the method is to identify both sides of equations (^ and ^ as Laurent 
polynomials, and to produce as many specializations of the variables that verify the equalities, as needed 
to imply these equations in full generality. 
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Zjj'p(27V; xi,. . . , xn-i,xn, ■ ■ ■ , X2N-i,x, y) 
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X2- 
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XN X2N-1 




= Ziij{2N + l;xi,. . . ,xn,xn+i, ■ ■ ■ ,X2N,x,y) 



XN+1 X2Ny 

Fig. 3: Partition functions for HTASM's 



4. 1 Laurent polynomials 

Since the weight of any vertex is a Laurent polynomial in the variables Xi, x and y, the partition functions 
are Laurent polynomials in these variables. Moreover they are centered Laurent polynomials, i.e. their 
lowest degree is the negative of their highest degree (called the half-width of the polynomial). In order 
to divide by two the number of non-zero coefficients (hence the number of required specializations) in x, 
we shall deal with Laurent polynomials of given parity in this variable. To do so, we group together the 
states with a given orientation (indicated as subscripts in the following notations) at the edge where the 
parameters x and y meet. 

So let us consider the partition functions ZQrp(4A^; X2n-i,x, y) and ZQy{AN] X2N-1, x, y), respec- 
tively odd and even parts of Zq^{AN] X2N-1, x, y) in x; ZQy{AN + 2; X2N, x, y) and ZQy{AN + 
2; X2N-, X, y), respectively odd and even parts of Zqj{AN + 2; X2n,x, y)\n x; Zpjrp(2Af; X2n-i,x, y) 
and Z^rp(2Af; X2n-ItX^ y), respectively parts with the parity of N and of A^— 1 of Z^j{2N; X2n~i, x, y) 



in x; and Z^rp(2A^ + 1; X2Ni x, y) and Zj|rp(2iV + \\X2Ni x, y), respectively parts with the parity of 
N - 1 andof iVof l;^2Ar,a;,y) in x. 

With these notations, the equations (H) and (^ are equivalent to the following: 

a{a)z'^r^{AN-X2N-i.x,y) = Z^r^{2N-X 2N-i,x, y)Z{N; X2n-i,x)Z{N; X2N-i,y), (10) 
a{a)z'^j{AN;X 2N-1, X, y) = z'^^{2N-X2N-i:X,y)Z{N-X2N-ux)Z{N-X2N-uy), (11) 
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Fig. 4: Partition functions for (q)QTASM of even size 



Quarter-turn Symmetric Alternating Sign Matrices 



55 



N, X, y)Z{N + 1;X2N, X, y)Z{N; X2n), (12) 



^HT 

7*1 



aia)Zf^j{4N + 2;X2N,x,y) = Z^j{2N+l;X2 N, X, y)Z{N + 1;X2N, x, y)Z{N; X2n)- (13) 

Lemma 7 Both left-hand side and right-hand side of equations f^^^^ are centered Laurent polynomials 
in the variable x, odd or even, of respective half-widths 2N — 1, 2N — 2, 2N, and 2N — 1. Thus to prove 
each of these identities we have to exhibit specializations of x for which the equality is true, and in number 
strictly exceeding the half-width. 

Proof: To compute the half-width of these partition functions, just count the number of vertices in the 
ice models, and take note that non-zero entries of the ASM {i.e. the first two orientations of Figure |l]) 
give constant weight (T(a^). Also, a line whose orientation changes between endpoints must have an odd 
(hence non-zero) number of these ±1 entries. □ 

4.2 Symmetries 

To produce many specializations from one, we shall use symmetry properties of the partition functions. 
The crucial tool to prove this is the Yang-Baxter equation that we recall below. 

Lemma 8 [Yang-Baxter equation] If xyz = a, then 



(14) 





z X 

The following lemma gives a (now classical) example of use of the Yang-Baxter equation. 
Lemma 9 



y- 

X- 



X- 

y- 



(15) 



Proof: We multiply the left-hand side by a(az), with z = axy. We get 



a{az) ^~ 



■:xz 
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□ 



The same method, together with the easy transformation 

5>Q ={a{az)+a{a^)) 

gives the following lemma. 
Lemma 10 



(16) 



y- 

X- 

y- 

X- 

y- 

X- 



cr(a^) + a{xy) x- 
a{a'^yx) y_ 

crjxy) X- 
aia^yx) y_ 

e^{xy) X- 



a{a?yx) y_ 



q-(a^) X- 

a{a?yx) y_ 

q-(a^) X- 

a{a?yx) y_ 



(17) 
(18) 
(19) 



We use Lemmas ^ and |l^ to obtain symmetry properties of the partition functions, that we summarize 
below, where m denotes either 2N or 2N + \. 

Lemma 11 The functions Z(N; X2n) ond Zj^^(2A^ + l]X2N,x,y) are symmetric separately in the 
two sets of variables {a;^, i < N} and {x^, i > A'^ + 1}, the function Z-^^{2N; X2N-1, a;, y) is sym- 
metric separately in the two sets of variables {xi, i < — 1} and {xi, i > N}, and the functions 
ZQ'p(2m; Xjv-i, x, y) are symmetric in their variables Xi. 

Moreover, Zqy{4N + 2; . . .) is symmetric in its variables x and y, and we have a pseudo-symmetry 
for ZQx(4iV; . . .) and Zyij{2N; . . .): 



ZQj{4:N;X2N~i,x,y) = 
Znj{2N;X2N-i,x,y) = 



a{a^)+a{xy) 

(j{a?yx) 
a{a?) + a{xy) 



'qt(4 



2N-1 



,y,x), 



cr(a^yx) 



— Znj{2N; X2N-i,y, x). 



(20) 
(21) 



Proof: For Z{N; . . .) and Zfj'p(m; . . .), the symmetry in two "consecutive" variables Xi and Xi+i is a 
direct consequence of Lemma ^ For ZQj(2m; . . .), we again apply Lemma ^ together with the easy 
observations: 

= and = _. ^ (22) 

which allow us to bring the Yang-Baxter triangle through the divalent vertices of Figure |[ 

For the (pseudo-)symmetries in (x, y), let us deal with ZQ^i^^'i ■ ■ ■)' other cases being similar or 
simpler. We use equation (^) to put together the lines of parameter x and y: 
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y X 



and then apply Lemma |T^. □ 

It should be clear that we have analogous properties for the even and odd parts of the partition functions. 
The next (and last) symmetry property, proved by Stroganov (H), appears when the parameter a equals the 
special value ajg = cxp(i7r/3). An elementary proof of this result has recently been given in 

Lemma 12 When a — ujq, the partition function Z{N; X2n) is symmetric in all its variables. 

4.3 Specializations, recurrences 

The aim of this section is to give the value of the partition functions in some specializations of the variable 
X or y. The first result is due to Kuperberg, the others are very similar. 

Lemma 13 [speciaUzation of Z; Kuperberg] If we denote 

A{xn+i,X2n\{xi,xn+i}) = W (j{axkXN+i) Yl (7{a'^XN+iXk), 

2<k<N N+l<k<2N 

A{xn+i,X2n\{xi,xn+i}) = Yl o-{axN+iXk) Y\. o-{a'^XkXN+i), 

2<k<N N+l<k<2N 

then we have: 

Z{N;ax.TSi+i,X2N\xi) = A{xn+i,X2n\{xi,xn+i})Z{N - 1\X2n\{xi,xn+i}), (23) 
Z{N;scxiTSi+:i_,X2N\xi) = A{xn+i, X2n\{xi,xn+i})Z{N - 1; X2n\{xi,xn+i})- (24) 
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Proof: We recall the method to prove equation (^3]). We observe that when xi — clxn+i, the parameter 
of the vertex at the crossing of the two lines of parameter xi and x m+i is a. Thus the weight of this vertex 
is a{aa) = cr(l) = unless the orientation of this vertex is the second on Figure [I] But this orientation 
implies the orientation of all vertices in the row xi and in the column xjv+i, as shown on Figure || The 
non-fixed part gives the partition function Z in size iV — 1, without parameters xi and xn+i, and the 
weights of the fixed part gives the factor A{. . .). 



XN- 



X2- 

Xi = axN+1- 



x\ — ax\- 

XN- 



X2- 



XN+1 



X2N 



XN+1 



X2N 



Fig. 5: Fixed edges for ( [23[ ) on the left and ( |24| ) on the right 
The case of (|4|) is similar, after using Lemma ^ to put the line a;Ar+i at the top of the grid. 



□ 



We will need the following application of the Yang-Baxter equation, which allows, under certain con- 
dition, a line with a change of parameter to go through a grid. 

Lemma 14 
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Proof: We iterative ly apply Lemma ^on the rows, and row by row 
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cx\. 



cx\. 
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—X 



□ 



Lemma 15 [specialization of Zy^^} If we denote 

A\j{xi,X2n\xi) = W (7{a^xiXk) W cj{axkXi), 

l<k<N N+l<k<2N 

Ah{xi,X2n\xi) = Y\. o-(a^a;feSi) Y\_ cr{axiXk), 

l<k<N N+l<k<2N 



A°jj{xn,X2n-i\xn) = Y\. cF{axkXN) W a{a'^XNXk), 

l<k<N-l N<k<2N-l 

'A°h{xn,X2n-i\xn) = Yl ^{axNXk) a{a^XkXN), 



l<k<N-l N<k<2N-l 

then for -k = ~i, *],*), and □ = 3, 5, *1 respectively, we have 

Zfjrp(27V+l;X2jv,a:,axi) = A]j{xi,X2N\xi)Z^r^{2N-X2N\xi,xi,x), (26) 
Zgrj,(2iV+ l;X2Ar,2;,axi) = A^H{xi,X2N\xi)Zl^r^{2N-X2N\xux,xi), (27) 
Zjjj(27V;X2Ar_i,a;,axN) = (T{axxN)A%{xN , X2n-i\xn)Z^j{2N - I; X2n-i\xn ,x,xn), (28) 
Z^rj,(27V;X2Ar_i,axN,j/) = a{axNy)A^H{xN ,X2n-i\xn)Zy[^{'2,N - 1; X2n-i\xn .U^xn)- (29) 

Proof: 

The proof is similar to the previous one, with the difference that before looking at fixed edges, we need 
to multiply the partition function by a given factor; we interpret this operation by a modification of the 
graph of the ice model, and apply Lemma It turns out that in each case, the additional factors are 
exactly cancelled by the weights of fixed vertices. 

To prove (^6|), we multiply the left-hand side by 



W aia^Xky), 



N+l<k<2N 



which is equivalent to adding to the line of parameter y a new line ay just below the grid; Lemma [14 
transforms the graph of Figure ^(a) into the graph of Figure ||(b). When we put y — axi, we get the 
indicated fixed edges, which gives as partition function 

Y[ <7^{aXkXi) Y\_ '7ia'^XiXk)Ziij{2N;X2N\xi,Xi,x). 

N+l<k<2N l<k<N 

Since a'^x^y = axkXi, the equation simplifies. To conclude, we observe that if we start with an edge 
going out from the crossing x/x2n (function ^jj-p) we get at the end the same orientation (function 
Z^j). □ 
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a^AT+l X2N 

(a) 

Fig. 6: Proof of 
Lemma 16 [specialization of ^Q'p] If we denote 

Aq{xi,X^^-y\x{) 

AQ{xi\Xm~l\x\) 




]^ a{(?XkXx)a(axiXk), 

\<k<m-\ 



J]^ a{(?xiXu)<y(axkXx), 

l<k<m.-l 

then for * = ^, ^, *1 and □ = *1, ^, ^, respectively, we have: 

ZQrj.(2TO;Xm_i,axi,?/) = a{axiy)AQ{xi,Xrn-i)ZQ^{27n-2]Xrn-i\xi,y,xi), (30) 
ZQrj,(2m;Xm_i,a;,axi) = a{axxi)AQ{xi; X„i^i\xi)Z^r^{2m - 2; Xra-i\xi,xi,x). (31) 

Proof: Similar to the proof of Lemma |5[ □ 

Remark 17 By using the (pseudo-)symmetry in [x, y), we may transform any specialization of the vari- 
able y into a specialization of the variable x. Moreover, by using Lemma ^^and (when a — loq) Lemma 
[72[ we obtain for Z, ZyYY '^'^d Zq^, 2N specializations. Now we have to compare them. 

4.4 Special value of the parameter a; conclusion 



When a = ujQ = exp(i7r/3), two new ingredients may be used. The first one is Lemma 12 , as mentioned 
in Remark [l^. The second one is that with this special value of a: 



(j{a) — iy{a?) a{a^x) = ^a{ax) = a{ax). 



(32) 
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which implies that the products appearing in Lemmas [l^, |l5| and [l^ may be written in a more compact 
way: 

A{xn+i,X2n\{xi,xn+i}) = o-{a) Y\_ cr{aXkXN+i), 
A{xn+i,X2n\{xi,xn+i}) = <j{a) Y\_ cr{axN+iXk), 



A]^{xi,X2n\xi) = Yl. cr{axkXi), 

l<k<2N 

Ah{xi,X2n\xi) = Y[ cr{axiXk), 



l<k<2N 



A%{xn,X2n~i\xn) = Yl <^{axkXN), 

l<fe<2JV-l 

'^h{xn,X2n^i\xn) = W a{axNXk), 



l<fe<2JV-l 



AQ{xi,Xrn-i\xi) = Y\_ a-'^{axiXk), 

l<fe<m-l 

AQ{xi,Xjn-i\xi) = Y\_ cr^iaXkXi). 



l<k<m-l 



Thus we get by comparing: 

A{xi,X2N\xi,x)AH{x^,X2N\xi) = a{axxi)AQ{x^,X2N\xi) 
A{xi,X2N\xi,x)A]j{xi,X2N\xi) = a{aXix)AQ{xi,X2N\xi), 

whence ( [lo| ) and (|ll]) imply that dl^ ) and (|l3) are true (in size 4iV + 2) for the 2N speciaUzations 
a; = a*^a;i (1 < i < A^). It is enough to prove ( 13 ) (Laurent polynomials of half-width 2N — 1), but we 
still need one specialization to get ( |l2| ) (half-width 2N). 

For (10) and ( 11 ), we observe the same kind of simplification 

A{xi,X2N-i\xi)(T{axxi)A%{xi,X2N-i\xi) = (T{axXi)AQ{xi, X2N-i\xi), 

whence ( [l^ and ( |l^ for the size AN — 2 imply that ( p^ and ( pi] ) are true for the N specializations 
X = axi, N < i < 2N — 1. We obtain in the same way the coincidence for the N specializations 
X — axi, N < i < 2N — 1. Thus we have 2N specialiations of x: it is enough both for (|l0|) (half-width 
2N - 1), and for (half-width 27V - 2). 
At this point, we have almost proved 

((^ and (|l]), in size 4iV) ((|l|) and (]^, in size AN + 2) ^ ((^ and (|llj), in size AN + A); 
almost, because we still need one specialization for (p^. 

We get this missing specialization, not directly for -^Q-p, ^qx' ^hx ^HT' ^'^^ original 
series Zqj(4A^ + 2; X2N, x, y) and ZYif{2N + 1; X2N, x, y): indeed if we set x — ay we may apply 
Lemma M . 
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Zqj{4:N + 2;X2N,ay:y) ^ a-{a) Yl '^iO'Xky)cria'^yxk)ZQj{4,N; X2n\x2N, 
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y 


















































































-ay 



XN+1 X2N 



Ziij{2N + l;X2N,ay,y) ^ \ (j{axky) G{a^yxk)\ Zy{-y:{2N;X2n\xn,xn-,xn) 

\l<k<N N+l<k<2N J 

This way, we get another point where is true, and thus, because we already have (|l3|), by difference 
we obtain that ( p^ holds for y — ax. 
This completes the proof of Theorem ^. 
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